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Abstract 



Two methodological troubles of the quantum theory of collisions are con- 
sidered. The first is the undesirable interference of the incident and scattered 
waves in the stationary approach to scattering. The second concerns the non- 
stationary approach to the theory of collisions of the type a + h ^ c + d. In 
order to calculate the cross section one uses the matrix element {cd\S\ah) of 
the S-matrix. The element is proportional to 5-function expressing the energy 
conservation. The corresponding probability |(c(i|S'|a6)p contains 5^ which is 
mathematically senseless. The known regular way to overcome the difficulty 
seems to be unsatisfactory. In this paper, both the troubles are resolved using 
wave packets of incident particles. 



1 INTRODUCTION 



An approach to the theory of scattering is known which may be called station- 
ary. The scattering is considered as a stationary process: there is a steady flux 
of particles incident on a potential V (the target). Scattered particles are also 
described by a steady flux. The state of the system is described by a vector 
which is constant in time. The vector is an eigenstate of the total Hamiltonian 
H = Hq + V which belongs to the continuous spectrum of H: Hipk = Ekipk- 
The eigenstate is known to be the superposition of the incident wave /^(x) and 
scattered wave Sf^{x): 



Here k = \/2mEk and x = (r, i^^ip), being the angle between k and x. The 
axis z is chosen to be parallel to the momentum k of the incident particle {x 
being its position). 

The stationary approach uses the probability flux or its density 



j{x) = — ^ 
2m 



V^/J*{x)^{x) - ^Ij*{x)V^{x) (2) 



instead of the usual probability |'?/'(a;)p, e.g., see [3], Ch. IV, § 29; [5J, Ch. II, 
§ 15, Ch. XI, § 95; [lOj, Ch. II.7. The probability flux generated by = I^+S^ 
(see Eq. (P) is 

Ik = — fv/i ■ Jr: - c.c.) + — (VS;- Sr- c.c. 



The cross section is determined by the ratio of the scattered flux (the second 
term in Eq. ([3])) to the incident one (the flrst term in Eq. ([3])), see Sect. |2]below. 
Besides these fluxes the total flux contains the interference terms (the third 
and fourth in Eq. ([3])). Their physical sense is unclear. It is unknown what 
contribution the interference terms may bring in the cross section. One may 
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conjecture that they must vanish if one replaces the plane incident wave 
by an incident wave packet, see p], Ch. X.5; [lOj, Ch. V, end of § 18. This 
conjecture is confirmed in Sect. [31 The wave packet is used which tends to the 
plane wave = exp{ikx) when packet dimension increases. Other ways of 
the packet introduction are possible. For example, Messiah [9] used a classical 
ensemble of small packets which have different impact parameters. Finally, 
averaging over the ensemble is carried out, see [9J, Ch. X. 

The stationary theory of scattering is nonrelativistic and is inapplicable, 
e.g., to the photon scattering (photon position operator and density of flux 
have no satisfactory definitions). The nonstationary approach is applicable 
to any process of the type a + b ^ c + d. It is based on the solution of the 
Schroedinger equation for the operator U{t, to) of evolution in time. Initially, 
at the moment to^ the system is in a state \E'(to) = {ab). At the moment t the 
system state is described by the vector \E'(t) = U(t,to)'^(tQ). The probability 
to find a final state \cd) at the moment t is equal to \{cd\U{t,to)\ab)\'^. The 
limit to — oo, t —>■ +00, i.e. S-matrix, is usually considered. 

The approach has the following trouble. The matrix element {cd\S\ab) is 
known to be proportional to the 5-function which expresses the total energy 
conservation: the total initial energy Ea+Eb is equal to the final energy E^+Ed 

{cd\S\ab) ^ 5{E), E = {E, + Ed)-iEa + Eh). (4) 

The probability | (c(i|S'|a6) p is proportional to the square 6"^ of this (5-function. 
This quantity does not exist mathematically, see [1]. Physicists gave to 6"^ an 
interpretation, see the end of Sect. IH but it cannot be recognized as satisfac- 
tory. Another resolution of this trouble is known, e.g., see [7], Ch. 14; [TU], 
Ch. VIII. It is presented in Sect. H] using the packet description of the initial 
state \ab). 

So a packet description of the incident particle allows us to resolve two 
troubles of the collision theory stated above. 
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2 Definitions of cross section 



The density of the incident flux F is deflned as the number of incident particles 
crossing per unit time a unit surface placed perpendicular to the direction of 
propagation. Let p be the number of particles per unit volume and v be the 
velocity of the incident particles. Then F = pv. If there is one particle in 1 
cm^, then F = v. 

Let us assume that the coordinate origin is in the center of a target. Let jr 
be the density of the probability flux of the scattered particles. The probability 
(or the number of particles) going through the area element rdipr smdd'd = 
r'^dQ during one unit of time is equal to jj.r'^dQ. This is the probability AA^ 
to detect the particle in the solid angle AQ during one unit of time (e.g., one 
second) 

AN = jyAVL. (5) 

The quantity AA^ may be related to the probability AW{t) of the particle 
detection in the solid angle at the moment t. One may assume 

AA^ = AW{t + 1 sec) - /\W{t) ^ ^AW{t) ■ 1 sec. (6) 

Usually another relation is assumed 

AA^ = W{t)/t. (7) 

Relations (El) and ([7]) coincide if time derivative d/SW {t) / dt is constant. One 
has the relation 

jyAVt = AN = AW{t) /t. (8) 

Therefore, the deflnition 

Aa = ^y^dn (9) 
of the cross section (see [3], Ch. XIII; [5], Ch. XI, § 95) is equivalent to 

Act = AN/F, (10) 

cf. |9], Ch. 10. 



3 Packets in stationary approach 

There is the conjecture that the description of an incident particle by a packet 
(instead of the plane wave) will turn into zero interference terms in the total 
flux (the third and fourth terms in Eq. ([3])). Let us give a conflrmation of this 
conjecture. 

The introduction of a packet implies that scattering is no longer a stationary 
process. Even if the scattering potential is absent, the (free) packet shifts and 
spreads. 

One possible setting of the problem of a packet scattering will be considered 
here (for another approach, see, e.g., [9], Ch. X). The definition of the cross 
section as the ratio of fluxes (see Eq. ([9])) will be retained though fluxes will 
not be stationary. The natural requirement is assumed (and ensured): in the 
limit when the packet turns into a plane wave the result should go to the usual 
stationary one. 

Consider a superposition 



of the H eigenfunctions H'^j: = Ek^j:. For \E'^ see ([1]), for see below. The 
superposition is not H eigenfunction, but the vector 

vi/(x,t)= j d'ke-^'^^'^^.ix)!^, (12) 

is a solution of the equation d'^{t)/dt = H'^it). 
The vector \l'(x, t) consists of two parts 

^(f,t) = I{x,t) + S{x,t), (13) 
I{x,t) = j rf=^A;e-*^'=*4y^^, (14) 

Six,t) = J d!^ke-''^^'i^A{^,^)e'^'/r. (15) 

The vector J(x, t) is the known description of the moving free packet (it is 
assumed that the spectrum oi H = Hq + V the same as the Hq spectrum). 



In order to calculate S{x, t) and the cross section, the following program is 
accepted. 

(a) The initial wave packet I{x, 0) will be chosen. 

(6) This determines the coefficients in Eqs. ( TTTl) . (IT^ and, therefore, S{x,t) 
may be calculated. 

(c) Absence of the interference terms will be verified. 

(d) Incident and scattering fluxes may then be found as well as the cross 
section. 

(a) Choice of /(x, 0). 

Consider the auxiliary wave function f{x) which is concentrated in a ball 
of the radius Rj (V/ being the ball volume). The ball will be named "support 



Note. In other words "support" is defined here as the volume outside which 
the function practically vanishes (or is unobservably small). In mathematics 
a different definition is accepted: the support is the volume outside which the 
function is exactly zero. 

If f{x) is spherically symmetric, then the average position J d?xxf*{x)f{x) 
is zero. Fourier transform of f{x) is also spherically symmetric and, therefore, 
the average momentum also equals zero. The average position of the shifted 
function f{x—a) = Id{x) is equal to a, average momentum being zero as before. 
One may verify that the function e^^^f{x) = I^{x) has average momentum p. 
At last, consider the shifted function Ip{x), i.e., the function 



Its average position is a and average momentum is p (the factor exp{—ipa) 



of fix) 



/p^(x-a)=e^P>^-'^)/( 



x-a)= 



(16) 



may be omitted). Let IjM.{k) be the Fourier transform of Ips{x). 



Let us assume 




0, outside the ball Vj 



1, inside Vj 



(17) 
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This means that /p'(x) = f{x) exp(zpx) is equal to exp(ip2f) inside Vj. When 
R] ^ oo the vector I^{x) tends to the plane wave whose wave function is equal 
to exp(ipx) everywhere. We have J d^x\Ij^{x)\'^ = Vf. there is one particle in 
the unit volume. 

The initial (at the moment t = 0) wave function of the incident packet is 
chosen to be equal to 

Ipa{S)=e'^^f{x-a). (18) 

The packet center a is placed at the point {—Ri) on the z axis which is parallel 
to p and passes through the potential center. This choice means that the 
interaction of the packet with the potential begins at the moment t = and 
stops at the moment T = 2Rj/v, v = P/mo (the packet dimension Rj is 
assumed to be much greater than the dimension i?^, of the volume in which 
the potential is concentrated). Till the moment t = and after the moment T 
the potential does not act on the packet and it is free. 

(6) Choosing I^{x) one may calculate S{x,t) (see f|T5|) ) using the Fourier 
transform f{p—k) exp{—ika) of /ps(x). If the packet I^{x) has a macroscopical 
dimension (e.g., 1 cm), then f{p—k) has a sharp maximum dXk = p. Therefore, 
Ek = Ep and \k\ = \p\. Let us calculate S{x,t) approximately using the 
expansions of Ej. and \k\ in Taylor series about the point p: 

Ek = Ep + {k-p)v, \k\ = \p\ + {k - p)u, (19) 
V = p/m, u = p/p, v = uv. (20) 

The scattering amplitude A{{}k, <^k) is simply replaced by A{{}p, ipp), -dk and 'dp 
being the angles between x and k, p, respectively, cf. Eq. ([1]). One obtains 

S{x, t) = A{^p, ^p)^expi{pr - Ept) j £ke'^'^-t^^"')e-''''^f{p- k). (21) 
After the change k' = p — k oi variables in j d^k . . . one gets 

S{x, t) = expi{pr - Ept)e"^ J d^k' f{k')e'^^'^\ 
b = —u{r — vt) + d, a = —uRj. 

6 



Here J cPk' f {k') exp{ik'b) is Fourier representation of /(&). Therefore 

S{x, t) = expi{pr - Ept)e-'^ f {-u[{r - vt) + Ri]). (22) 

Remind that f{h) = 1 when \b\ < Rj, i.e. | — u[{r — ft) + Rj]\ < Rj or 

\r -vt + Ri\ < Ri, (23) 

because \u\ = 1 and / is spherically symmetric. The support of S{x,t) is 
determined by inequality (!23|) . Let us discuss it. 

One may assume that r ^ Rj: the packet dimension Rj is much less than 
the distance between the target and the detector. Then fl25]) may be replaced 
by a simpler inequality 

|r - vt\ < Ri. (24) 

It determines the movement of the scattered packet. Let us compare it with 
the movement of the initial incident packet Ij;{x) = f{x) expipx, see Eq. ( JTSl) . 
The latter is the shift x ^ x — vt along the z axis, if the spreading is neglected, 
(e.g., see [I2], Ch. 10.4): 

I^{x,t) = I^{x — vt,0) = f{x — vt) expip{x — vt), v=p/m. (25) 

This means the shift of the packet support, i.e., the ball of the radius Rj 
(the phase factor expipvt is inessential). The shifted ball is described by the 
inequality 

\x-vt\<Ri. (26) 

The shifted ball is not a spherically symmetric region (excluding the case 
when the packet centre coincides with the coordinate origin, i.e., potential 
centre). Meanwhile, inequality (I24p describes a spherically symmetric region 
at all times: flM|) does not contain angles i), if of the vector x but contains 
only \x\ = r. At fixed t the support region is the spherical layer between the 
spheres of the radii vt — Rj and vt + Rj. The thickness of the layer is equal to 
2Rj. As if: increases, this layer "inflates" preserving its thickness. 
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When t is small and \r — vt\ = |r|, then inequality fl23l) is not fulfilled 
and S{x,t) = 0: the scattered wave appears in the detector only when t is 
sufficiently large. 

(c) We obtain that the scattered packet is in the spherical layer described 
above, while the support of the incident packet moves along the z axis. The 
moving packets' supports do not practically intersect if r is sufficiently large. 
Therefore, S{x,t)I^{x,t) =0: is zero where / is nonzero and vice versa. 
So the interference of the incident and scattered waves is absent (with the 
exception of small values of the angle ■&). 

(d) The used expressions f l2^ and (1221) for incident and scattered wave 
packets differ from the corresponding and Sk waves in Eq. ([T]) only in one 
respect: the former have an additional factor / which is equal to unity inside 
the moving packets and vanishes outside them. Therefore, the incident and 
scattered ffuxes inside packets are the same as in the stationary case. However, 
these ffuxes are nonstationary: their supports move in space. One may retain 
the previous deffnition iQ for the cross section having in mind that the ffuxes 
F and jy in ([9]) are intrapacket ones. 

Note once more that the used packet's description of scattering turns into 
the ordinary stationary one when Rj oo (the condition r ^ Rj being 
implied) . 

4 Packets in nonstationary approach 

The nonstationary approach to the collision theory uses the evolution operator 
U{t, to) (interaction or Dirac picture is in mind). For the reaction a + 6 — > c+d 
one must calculate the matrix element of the type {cd\U{t,to)\ab) . Let us 
assume that the initial state \i) = \ab) is the product of packets |a), \b), see 
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Sect. 131 For example 

\a) = J d^k\k)a{k). 

The packets have finite supports in the coordinate space. The supports are sup- 
posed to be of macroscopically large dimensions and, therefore, packet spread- 
ing may be neglected, see [6], Ch. 3.1; [12], Ch. 10.4. 

The particles do not interact if their supports do not intersect. So the 
interaction lasts during a finite interval T of time. In the following I let to = 
-T/2, t = T/2. 

Let us consider the matrix element of the expansion of U{T/2, —T/2) in 
the perturbation series 

rT/2 

Uf,{T) ^ {f\U{T/2, -T/2)\^) = + z / dte^^^f-^^^^mtjz) 

J -T/2 

pT/2 nti 

/ citie'(^/-^-)*^ / dt.S^--^^'^'^ {mtJm) {m\HtJi) +(27.) 

Cf. [S], Chs. 1.2 and 1.3; [11], Ch. 11.6. Here Ef denotes the final total energy 
Ef = Ec + Ed- Analogously, Ei = Ea + Eh, {Ea and Eh being average energies); 
Em is the total energy of the intermediate state |m); H-^^ is the interaction 
Hamiltonian in the Schroedinger picture. One may suppose that {f\i) = 0. 
The second term in Eq. f l27l) is proportional to 

27T6TiE) = dte""' = '^^ELll^ E = Ef- E,. (28) 

J~T/2 E 

The third term contains the integral over ^2 

J -T/2 

The contribution to (1291) originating from the lower limit —T/2 tends to zero as 
T ^ oo due to fast oscillations of exp i{Em — E^T /2. For a strict proof of this 
statement one must use the packet description of \i) and Riemann-Lebesque 
lemma, see, e.g., [H], Ch. 11 after Eq. (11.165). Neglecting this contribution 
one obtains that the remaining integral over ti is equal to 2n6T{E), Eq. (!28l) . 
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Analogously, one may argue that in all orders of the expansion (127|) U fi (T) is 
proportional to SxiEf — Ei) 

Uf,{T) = 6T{Ej-Ei){f\R\z), (30) 

where {f\R\i) ceases to depend upon T as T ^ oo. Note that St{E) S{E) 
as T — s> oo. In this limit Eq. (1271) turns into 

Sf, = 2-KS{Ef-E,){f\R\t), (31) 

where S is the S'-matrix. One gets that the probability |(/|f/(oo, — oo)|2)p is 
proportional to 6"^. The square of the 5-function has no mathematical sense, 
see [1], part 111, Sect. 12.5. 

However, for the cross section we need probability in unit of time, see 
Sect. H It may be defined as |f//,(T)|VT. It follows from Eqs. ([30]) and (|28|) 
that 

l/-/._.m|2/T ^ 1^ 

TE^ 



U,m\VT.'J^. (32) 



The r.h.s. of Eq. ([321) tends to 2-k5{E) as T ^ oo, not to 5'^{E), see [7], Ch. 2, 
Eq. (8.19). So we obtain the following value for the probability in unit time: 

hm \Uf.{T)\^/T = 27i5{Ef - Em\m^- (33) 

1 — *oo 

The probability in unit time may be defined in a different way, namely as 
d\Uf,{T)\'^/dT. We have 

d .^M9 d 4sin2ET/2 2 , 
— \Ufi{T)\^ ~ — ' = - sin ET 2. 34 

In the limit T — > oo one gets in the r.h.s. of Eq. (1331) the 5{E) function as 
above, cf. m], Ch. 11, Eq. (11.91). 

The usual way to calculate the cross section is to start with the S-matrix 
element (/|f/(oo, ~oo)\i). The reason is that it is relativistic invariant, allows 
renormalization etc., unlike (/|f/(t,to)K)- But the probability |(/|>S'|z)p is pro- 
portional to 6'^{Ef — Ei) and this is senseless. The trouble is usually overcome 
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in the following manner (see, e.g., [8j, Ch. 1.2; [llj, Ch. 14.1; [2], Ch. 5, § 37; 
^). In the product 6{E)6{E) one of the 5-functions is replaced by 6{0) because 
of the presence of the other 5-function. Basing on the representation 



(5(0) is replaced by T/27T. In order to obtain the probability in unit of time 
one divides K/I^l^)!^ by T. So one gets 



which is a sensible result. 

However, the argumentation is not satisfactory: S{E) does not depend on 
T unlike St{E), see Eq. ([28]) {S{E) is the hmit of (EHD as T ^ oo). Instead of 
(!35|) one may use the representation 



and analogously obtain 5(0) = 2T/27r instead of 5(0) = T /2ti. 

Nevertheless, the resulting "probability in unit time" obtained in the books 
cited above coincides with the r.h.s. of Eq. fl33p obtained here. 



As has been expected, the interference of the incident and scattered waves 
in the stationary theory of scattering is absent if the waves are described by 
packets. The cause is the nonintersection of the packets' supports (excepting 
the limitingly small scattering angles). 

Consideration of the 5^-trouble arising in the nonstationary approach needs 
the determination of the "probability in unit time" Wfi{T). In Sect. HI for this 
purpose the evolution operator U{T/2, — T/2) was used, Wfi{T) being defined 
as 




(35) 



\{f\S\^)\yT ^ 6{E) 




5 CONCLUSION 



Wf,{T) = \{f\UiT/2,-T/2)\t)\yT. 



(36) 
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In the limit T — > oo one gets (see Sect. Hj), 

\imWfi{T) 6{E), E = Ej - Ei. (37) 

The result is free of the 5"^ trouble. 

Another way to get "probability in unit time" is given in books on quantum 
field theory, e.g., see [S], [H], [2], \M- first, one considers the limf/(T/2, —T/2) 
as T ^ oo, i.e., the S'-matrix (in our way the limit T ^ oo is carried out later, 
see (l36l) . (1371) ). The S'-matrix elements {f\S\i) are proportional to S{E). The 
corresponding probability is proportional to S'^{E), which is sense- 

less. In the books 5^ is treated in an unsatisfactory manner presented and 
criticized at the end of Sect. HI 

Both the ways give the same result, Eq. fl33l) . Here the satisfactory way of 
getting the result is considered. 

References 

[1] p. Antosik, J. Mikusinski, R. Sikorski. Theory of Distributions, part III, 
Sect. 12.5. Elsevier Scientific Publishing Company, Amsterdam, 1973 

[2] A. Baldin, V. Goldanskii, I. Rozenthal. Kinematics of Nuclear Reactions. 
Pergamon Press, Oxford, 1961. 

[3] D. Blokhintsev. Quantum Mechanics. Reidel Publ. Co., Dordrecht- 
Holland, 1964. 

[4] N. Bogoliubov, D. Shirkov. An Introduction to the Theory of Quantized 
Fields, Ch. IV, § 24.5. John Wiley, New York, 1980. N.N. Bogoliubov, D. 
Shirkov. Quantum Fields, Ch. V, § 21.1. Benjamin / Cummings, London, 
1983. 

[5] A. Davydov. Quantum Mechanics. Pergamon, New York, 1976. 



12 



[6] M. Goldberger, K. Watson. Collision Theory. John Wiley, New York, 
1964. 

[7] W. Heitler. The Quantum Theory of Radiation. Clarendon Press, Oxford, 
1954. 

[8] P. Matthews. The Relativistic Quantum Theory of Elementary Particle 
Interactions. Rochester, New York, 1957. 

[9] A. Messiah. Quantum Mechanics, V. 1. North Holland, Amsterdam, 1961. 

[10] L. Schiff. Quantum Mechanics. McGraw-Hill Book Company, New York, 
1955. 

[11] S. Schweber. An Introduction to Relativistic Quantum Field Theory. Row, 
Peterson, New York, 1961. 

[12] J. Taylor. Scattering Theory. The Quantum Theory of Nonrelativistic Col- 
lisions. John Wiley, New York, 1972 



13 



